Photon localization and Dicke superradiance in atomic gases 
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Photon propagation in a gas of N atoms is studied using an effective Hamiltonian describing 
photon mediated atomic dipolar interactions. The density P(T) of photon escape rates is determined 
from the spectrum of the N x N random matrix I\j = sin(xij)/xij, where Xij is the dimensionless 
random distance between any two atoms. Varying disorder and system size, a scaling behavior is 
observed for the escape rates. It is explained using microscopic calculations and a stochastic model 
which emphasizes the role of cooperative effects in photon localization and provides an interesting 
relation with statistical properties of "small world networks". 
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We address the question of photon localization result- 
ing both from disorder and cooperative effects (multi- 
atomic coherent emission) in three-dimensional cold 
atomic gases. This localization shows up as an overall 
decrease of photon escape rates and our purpose is to 
investigate the roles played by disorder and by coopera- 
tive effects. For weak disorder, incoherent spontaneous 
emission by independent atoms is expected. For stronger 
disorder, cooperative effects become important and lead 
to vanishing escape rates, i.e. to photons trapped in the 
gas for very long times. We show that localization occurs 
as a smooth crossover between these two limits and not as 
a phase transition like for Anderson localization [1]. We 
propose a stochastic description of the photon emission 
process, which agrees quantitatively with our numerical 
results and explains the nature of the crossover. This 
leads to the conclusion that photon localization in atomic 
gases is primarily determined by cooperative effects and 
not by disorder. 

We consider a collection of N identical atoms at rest, 
taken to be degenerate two-level systems respectively de- 



noted for the atom i, by \gi) = \j g 



0,m o 



0) and 



\ e i) = lie = l,m e = 0, ±1} for the ground and excited 
states, j is the total angular momentum and m is its 
projection on a quantization axis. The energy separation 
between the two levels, including the radiative shift, is 
Hluq and the natural width of the excited level is HTq. 
Atoms randomly placed at positions r^, are coupled to 
the electromagnetic field E through their dipole opera- 
tor dj. The corresponding Hamiltonian is, 



N 



N 



i=l 



ke 



where a ke is the creation operator of a photon of wave 
vector k, LOk — c|k and polarization e. We assume that 
only one photon is present. For resonant scattering, trac- 
ing over the photon degrees of freedom leads to the ef- 



fective atomic Hamiltonian 
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where = |ej)(<7j| is the atomic raising operator, 

S7 = {S^Y and S z = J2i=i S zi witn S zi = |ei)(e»|. 
The potential Vy = /3y — vjij is a random and complex 
valued quantity, specified by 
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where ko — ujo/c and n 



r,-| is the distance be- 



tween any two atoms. The quantities p and q depend 
on the atomic transition. For Am = m e — m„ = 0, 



po = sin tfij 
P± = 1(1' 

Oij - 



and qo = 1 — Scos 2 ^ and for Am = ±1, 
I- cos 2 9ij),q± = -|(3cos 2 #ij — 1), where 
= cos _1 (£ • iij) and is a unit vector defined along 
the direction joining the two atoms. Expressions (2) and 
(3) are well known [2] and correspond to an instanta- 
neous photon exchange between an initially excited atom 
i coupled to the j-th atom. This description holds for 
distances between atoms that are small compared to the 
coherence length of the light emitted by a single atom 
[3]- 

The possible escape rates T of a photon propagating 
in an atomic gas, are obtained from the time evolution 
of the ground state population (G\p\G), 

"' - ' ~ ^ -3ipSt\G) (4) 



dt (G\ P \G) = r J2^(G\S7pst\G) 



where \G) = \g±, <?2, • ■ ■ , 5at) and p(t) is the reduced 
atomic density operator [4] . The real and symmetric ma- 
trix 7jj can be diagonalized by an orthogonal transfor- 
mation. The eigenvalue equation 
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together with the collective operators = 
E*Li«i* ) Sj t . enable to rewrite (4) as (G\p\G) = 
T J2k=i r k(G\S^pS+\G). This last form allows to 
interpret the eigenvalues {T^} as the escape rates and 
the eigenf unctions {uj; } as photon modes providing 
the directivity of the angular emission [4]. It results 
from (4) that the IVs are independent of the effective 
dipole-dipole interaction /3y so that Van der Waals 
dephasing does not play a role [5]. 

The average density of escape rates P(T), normalized 
to unity, is defined by P(T) = -(1/tt) Imi?(z = r + iO+), 
with 
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The average 7TT is taken, at fixed density, over spatial 
configurations of the atoms. For Gaussian ensembles of 
random matrices [6], P(T) obeys a semi-circle law. Here, 
as we shall see, the behavior is very different. 

The random matrix 7y depends on the distances 
between atoms and on the angles Oij. We expect local- 
ization properties to depend on rather than on 9^. 
We therefore consider the scalar model [7] obtained from 
(3) by averaging 7^ over 6^ and thus given by the N x N 
random matrix [8], 



(7) 



where Xij — k^rij are interatomic distances expressed in 
units of the wavelength A = 2n/kQ. We have checked 
that the vectorial and scalar behaviors of P(r) are qual- 
itatively the same [9]. This point is important since one 
may think that the scalar case (7) is the far-field limit 
obtained by dropping in (3) the near-field terms respon- 
sible for cooperative effects. This is not the case and Ty 
is well defined for iy = so that TrfTjj] — N . The 
eigenvalues of T^ are nonnegative since the 3d Fourier 
transform of sinc|a;| is S(\k\ — 1) > 0. 

We now consider N atoms enclosed in a cubic volume 
L 3 , with a uniform density n. The disorder strength 
is defined by the dimensionless parameter W = l/k^l 
where I = 1/na is the elastic mean free path and a ~ A 2 
is the resonant scattering cross section [10]. Introducing 
the number N± = (fcoL) 2 /4 of transverse photon modes, 
leads to W = f^^j- 

Characteristic behaviors of P(T) for different values of 
W and size a = L/X are displayed in Fig.l. For a dilute 
gas (W <C 1) we recover the single atom limit Tij — > Sij, 
namely P(T) is narrowly peaked around T = 1 (in units 
of To) as expected from resonant scattering of a photon 
by a single atom (Fig.l. a). For stronger disorder, P(T) 
becomes broader and shifted towards lower values of T 
(Fig.l.b). Eventually for large enough disorder, most of 
the eigenvalues get close to T = (Fig.l.c). Such a van- 
ishing escape rate corresponds to photons localized in the 
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FIG. 1: Behavior of P(T) for different values of disorder W 
and size a — L/X, with N = 216. (a) low disorder (b), (c) 
larger disorder and (d) Dicke limit. 



atomic gas. By increasing further W, at fixed number N 
of atoms, yet another behavior shows up for x\j <C 1 
(Fig.l.d) where P(T), obtained from the escape matrix 
with all entries equal to one (Tij = 1), has two eigenval- 
ues. One at r = is the (N — l)-degenerate subradiant 
mode and the second V = N is the non-degenerate super- 
radiant mode. This is the Dicke limit [11] reached when 
atoms are enclosed in a volume much smaller than A 3 . 
Using (6), we obtain the density: 

P(T) = ^5(T) + ^S(T-N) . (8) 

A quantitative characterization of the behavior of 
P(r) is obtained by considering the relative number of 
states Jj 00 dT P(T) having an escape rate larger than 1. 
We then introduce the conveniently normalized function 
C(a, W) defined between and 1 by 



C(a,W) = 1 



dT P(T) 



(9) 



C(a, W) thus defined, measures the relative number of 
states having a vanishing escape rate. At finite size, we 
expect C(a, W) to have a scaling form [12], namely : 



d\nC(a,W) 
din a 
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(10) 



whose solution C(a,W) is a function of a/£(W) alone. 
We have verified this scaling behavior over a broad range 
of size and disorder. For a > 1, the results displayed in 
Fig. 2, collapse on a single curve (Fig. 3) when plotted as a 
function of the parameter 2waW = ir 2 N/N±. The Dicke 
limit (Fig.l.d) is reached for small volumes, L <C A, and 
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FIG. 4: Behavior of C(N) in the Dicke limit. 



FIG. 2: Behavior of C as a function of the size a = L/X > 1, 
for different disorder strengths W . 




FIG. 3: All points represented in Fig. 2 collapse on a single 
curve as expected from (10). The solid line is a fit of (14) 
and the dashed line is the asymptotic behavior obtained from 
(12). 



N± ~ 1, so that the scaling parameter becomes N/N± ~ 
N. Using (8) and (9), leads to C(a,W) = 1 - (2/N) 
displayed in Fig. 4. 

To explain this scaling behavior, we study first the 
limit N < Nj_ which can be understood from sim- 
ple considerations on the superradiant rate. The elec- 
tric held E created by N excited atoms within the vol- 
ume L 3 , is obtained from the electromagnetic energy 
NTiluq = L 2 (c/T)E 2 /8ir. The superradiant escape rate T 
may also be obtained by assuming that under the action 
of the field E, each atom performs half a Rabi oscillation 
over a time 1/T = h/Ed. The dipolar matrix element d 

is related to the spontaneous emission rate Tq = f ■ 
Altogether, T is specified by the parameter jj- as 
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E is delocalized over the atoms i.e. for N/N± < 1, 
emphasizes that the initial linear behavior of C(N/N±) 
(Fig. 3) is essentially determined by cooperative effects. 
In the opposite limit N 3> Nj_, the n-th order cumulant 
of P(r), is (l/JV)Trrjr = ^(N/N^Y 1 - 1 /{n + 2), and a 
resummation leads to the asymptotic behavior [9], 
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for r < N/N± and P(T) = otherwise, so that asymp- 
totically, C(N/N±) = 1 - 3(N±/N) (Fig.3). 

To interpolate between the two previous limits, we con- 
sider a mapping of the cooperative emission of randomly 
distributed atoms onto a stochastic Markov process on a 
graph. To define it, we start from the Dicke limit [11] 
whose escape rate matrix = 1, is the adjacency ma- 
trix of a complete graph having spectral density (8). In 
that limit, the Hamiltonian (2) commutes with the collec- 
tive spin operators 5' ± = J2iLi so that Dicke states 
1 5, m), eigenstates of S 2 = ±{S+S~ + S~S+) + \S 2 Z and 
S z , are eigenstates of H e . The emission of a photon, a 
process which keeps S unchanged and shifts m by one 
unit, can then be described as a one-dimensional Markov 
process, P/v(0), with equiprobable jumps (1/2) between 
neighboring m-states. Away from the Dicke limit, the ex- 
pressions of the collective spin operators for scalar waves, 



are obtained by taking 



Xk) 



(5), so that 



The random phases prevent from 



having the previous angular momentum algebra so that 
H e does not commute with Yet, we can still use the 
Dicke states basis, but now a photon emission is a process 
where both S and m change. This can be described as 
a modified Markov process where S'-changing events are 
accounted by adding random jumps to non-neighboring 
m-states with a probability e. The corresponding Markov 
process P/v(e) is defined by the N x N matrix 



Pjv(e) = (1 - Ne)P N (0) + e\e N )(e N \ 



(13) 



This expression which is correct when the electric field 



where Pjv(0)i,j±i = 1/2. The escape rate of a photon 
is set by the inverse mean hitting time 7/v(e) defined as 
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the average over all sites of the number of steps needed 
to reach a given assigned site [13]. The scaling function 
(9) can be written as C(e,N) = eNT%(e)/T N (Q) and for 
large N, C depends on s = (7/ tanh7) — 1 alone, namely 



where 7 = a/ eN 3 /2 [9, 13]. The parameter s is the num- 
ber of shortcuts between non neighboring m-states in- 
duced by emission processes which do not leave S un- 
changed. A plot of (14) is shown in Fig. 3 and compared 
to the numerical data. 

Previous expressions indicate that photon escape rates 
are primarily determined by cooperative effects and not 
by disorder. Moreover, (14) displays a smooth crossover 
between delocalized and localized photons and not a dis- 
order driven localization transition. To emphasize this 
point, we compare the two parameters, N/N± which 
specifies the strength of cooperative effects, and the di- 
mcnsionless conductance g = N]_/N which determines 
the strength of disorder [10]. These two quantities are 
plotted in Fig.5 for a fixed, but large N. For Nj_ > TV 2 / 3 , 
g is large so that disorder effects are small. In the oppo- 
site limit, we have N/N± ;§> g so that cooperative effects 
prevail disorder in a regime (g <C 1) where disorder may 
precisely be strong enough to lead to disorder driven lo- 
calization. 




FIG. 5: Comparison between cooperative effects determined by 
the parameter N/N± and disorder described by g = N±/N. 

To summarize, we have shown that escape rates of pho- 
tons propagating in a 3d atomic gas, are characterized by 
a scaling behavior determined by the parameter N/N±. 
For N/N± <C 1, cooperative effects are negligible and 
photons are emitted in spontaneous and incoherent pro- 
cesses. For larger values of this parameter, cooperative 
effects set in and the overall photon escape rate becomes 
very small. Eventually, for fixed N, we reach the Dicke 
limit. The crossover between these behaviors is smooth 
and a disorder induced localization phase transition is 
unlikely to take place. 



These results are well described by the modified 
Markov process (13). Surprisingly enough, this process 
also provides an accurate description and an interesting 
link [13] to the recently widely studied "small world net- 
works" [14], tuned to be intermediate between regular 
networks with long (oc N) chemical lengths and random 
networks with short (oc In N) ones. The existence of a 
crossover rather than a phase transition, between regu- 
lar and random networks has been thoroughly investi- 
gated [15]. Noting that "small world networks" appear 
to be well suited to achieve synchronization of non lin- 
ear oscillators [16] and that photon cooperative emission 
results from the synchronization of the atomic dipoles in- 
duced by long range atomic correlations, this connexion 
becomes even more interesting and relevant. Finally, the 
analysis we present in this letter may suggest a different 
approach and new protocoles for experiments on photon 
localization in cold atomic gases. 
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